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Abstract 

We consider the determination of the solar neutrino oscillation parameters Am^ 
and 9\2 by studying oscillations of reactor anti-neutrinos emitted by nuclear power 
plants (located mainly in France) with a detector installed in the Frejus underground 
laboratory. The performances of a water Cerenkov detector of 147 kt fiducial mass 
doped with 0.1% of Gadolinium (MEMPHYS-Gd) and of a 50 kt scale liquid scintillator 
detector (LENA) are compared. In both cases 3a uncertainties below 3% on Am^ and 
of about 20% on sin 2 #12 can be obtained after one year of data taking. The Gadolinium 
doped Super-Kamiokande detector (SK-Gd) in Japan can reach a similar precision if 
the SK/MEMPHYS fiducial mass ratio of 1 to 7 is compensated by a longer SK-Gd 
data taking time. Several years of reactor neutrino data collected by MEMPHYS- 
Gd or LENA would allow a determination of Am|i and sin 2 612 with uncertainties of 
approximately 1% and 10% at 3a, respectively. These accuracies are comparable to 
those that can be reached in the measurement of the atmospheric neutrino oscillation 
parameters Am^ and sin 2 623 in long-baseline superbeam experiments. 



a Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
Bulgaria. 
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1 Introduction 



The experiments with solar [1-3], atmospheric [4], reactor [5,6] and accelerator neutri- 
nos [8, 9] have provided during the last several years compelling evidence for existence of 
neutrino oscillations caused by nonzero neutrino masses and neutrino mixing. The data 
imply the presence of 3-u mixing in the weak charged lepton current (see, e.g., [10]): 

3 

vil = ^2 Ul i ^' L ' l = e ^^ r > (!) 

3=1 

where are the flavour neutrino fields, z/jl is the field of neutrino i/j having a mass mj and 
U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [11], U = £/pmns- All 
currently existing //-oscillation data, except the data of the LSND experiment [12], can be 
described perfectly well assuming 3-u mixing in vacuum and we will consider this possibility 
in what follows. 1 

The PMNS matrix can be parametrized by three angles and, depending on whether the 
massive neutrinos Uj are Dirac or Majorana particles, by one or three CP-violation (CPV) 
phases (see, e.g., [16]): 

t/pMNs = v(e 12 , e 13 , e 23 , 5) dia g (i, e'f, e *f ) (2) 

where V(9i2, #13, 623, 8) is a CKM-like matrix, 8 is the Dirac CP-violating phase and a, f3 are 
two Majorana CPV phases [17,18]. If we standardly identify Am| = Am^ = m 2 — m 2 > 0, 
where Amg drives the solar neutrino oscillations, then Am\ = Am^ = Am§ 2 , O23 = 9 A and 
9 12 = 00, Am^, 9a and 9 & being the z/-mass squared difference and mixing angles responsible 
respectively for atmospheric and solar neutrino oscillations, while #13 is the mixing angle 
constrained by the CHOOZ experiment [7]. The existing neutrino oscillation data allow us 
to determine Aro^, jAm^j, sin 2 9 12 and sin 2 2# 2 3 with a relatively good precision and to 
obtain rather stringent limits on sin 2 9 13 (see, e.g., [19-21]). The best fit values and the 
95% C.L. allowed ranges of Am^, sin 2 #i 2 , lAmjjjJ and sin 2 2^23 read: 

Am^ = 8.0 x 10~ 5 eV 2 , sin 2 21 = 0.31 , 

I Am 2 I = 2.2 x 10" 3 eV 2 , sin 2 2 9 23 = 1.0 ; 

(3) 

Am 2 x = (7.3 - 8.5) x 10~ 5 eV 2 , sin 2 9 12 = (0.26 - 0.36) , 
I Am 2 x I = (1.7-2.9) x 10~ 3 eV 2 , sin 2 2 9 23 > 0.90. 

A combined 3-v oscillation analysis of the global data gives [19,21] 

sin 2 13 < 0.027 (0.044) at 95% (99.73%) C.L. (4) 

Using the recently announced (but still unpublished) data from the MINOS experiment [9] 
in the analysis leads to a somewhat different best fit value and 95% allowed range of lAm^l 

1 In the LSND experiment indications for — ► v e oscillations with Ato^snd — 1 °V 2 were obtained. 
The minimal A-v mixing scheme which could incorporate the LSND indications for i/-oscillations is strongly 
disfavored by the data [13]. A ^-oscillation explanation of the LSND results might be possible assuming 5-v 
mixing [14]. The LSND results are being tested in the MiniBooNE experiment [15]. 
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and to a somewhat more stringent limit on sin 2 13 [21]: (Am 2 -,] = 2.6 x 10 3 eV 2 , (Am^-J = 
(2.2 - 3.0) x 10- 3 eV 2 , and sin 2 13 < 0.025 (0.040) at 95% (99.73%) C.L. 

In spite of the enormous progress made in establishing the existence of neutrino oscil- 
lations driven by non-zero neutrino masses and mixing and in determining the pattern of 
neutrino mixing and the values of the two neutrino mass squared differences, our knowledge 
and understanding of neutrino mixing is rather limited at present (see, e.g., [10] for a detailed 
discussion of the current status of our ignorance about neutrino mixing). Future progress 
in the studies of neutrino mixing requires, in particular, the knowledge of the precise values 
of the parameters which drive the solar and the dominant atmospheric neutrino oscillations, 
Amjj, sin 2 # 12 , Am^ and sin 2 #23 (see, e.g., [22]). The high precision measurement of these 
parameters is one of the main goals of the next generation of neutrino oscillation experi- 
ments. In the present article we discuss the possibility of a high precision measurement of 
the solar neutrino oscillation parameters Ara^ and sin 2 6*12 in an experiment studying the 
oscillations of reactor anti-neutrinos V e with a "large scale" detector located in the Frejus 
underground laboratory in France. 

The existing data allow a determination of Am 2 x and sin 2 9 12 at 3cr with an error of 
approximately 11% and 25%, respectively. These parameters can and will be measured 
with higher precision in the future. The data from phase-Ill of the SNO experiment [3] 
using 3 He proportional counters for the neutral current rate measurement could lead to a 
reduction of the error in sin 2 9 12 to 21% [23,24]. If instead of 766.3 t yr one uses simulated 
3 kt yr KamLAND data in the same global solar and reactor neutrino data analysis, the 3a 
errors in Am^ and sin 2 9 12 diminish to 7% and 18% [24]. The most precise measurement 
of Am^, discussed so far in the literature, could be achieved [23] using Super-Kamiokande 
doped with 0.1% of Gadolinium (SK-Gd) for detection of reactor v e [25]: the SK detector 
gets the same flux of reactor v e as KamLAND and after 3 years of data-taking, Am^ could 
be determined with an error of 3.5% at 3er [23]. A dedicated reactor v e experiment with 
a baseline L ~ 60 km, tuned to the minimum of the v e survival probability, could provide 
the most precise determination of sin 2 9\ 2 [26]: with statistics of ~ 60 GW kt yr and a 
systematic error of 2% (5%), sin 2 12 could be measured with an accuracy of 6% (9%) at 
3o- [24]. The inclusion of the uncertainty in 9 13 (sin 2 9 13 <0.05) in the analyzes increases the 
quoted errors by (1-3)% to approximately 9% (12%) [24]. The improved determination of 
Amjj and 9 12 with KamLAND or dedicated post-KamLAND reactor neutrino experiments 
has been studied previously also in Refs. [27-31], whereas the potential improvements of the 
precision on these parameters from future solar neutrino experiments has been investigated, 
e.g., in Refs. [24,26,32,33]. 

MEMPHYS (MEgaton Mass PHYSics) [34] is a project for a mega ton scale water 
Cerenkov detector located in the Frejus underground laboratory at the border of France and 
Italy. It is similar to the UNO [35] project in the US and the future Hyper-Kamiokande [36] 
detector in Japan. Such detectors allow for a broad range of physics studies like nucleon 
decay, long-baseline accelerator neutrino oscillations, super nova neutrino detection, and os- 
cillations of solar and atmospheric neutrinos. The MEMPHYS detector is considered as a 
far detector for neutrino beams produced at CERN located at a distance of 130 km from 
Frejus (see, e.g., Ref. [37]). A recent civil engineering pre-study indicates that MEMPHYS 
could be built with existing techniques as a modular detector consisting of three (up to five) 
modules (shafts), each having a fiducial mass of approximately 147 kt. 
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In the present paper we consider the possibility that the water in one module of MEM- 
PHYS is doped with 0.1% of Gadolinium (MEMPHYS-Gd), as it has been proposed originally 
for Super-K [25]. This allows a very efficient detection of electron anti- neutrinos through 
the reaction v e + p — > e + + n since the neutron can be tagged due to the high absorption 
cross section on Gadolinium. One module of MEMPHYS (147 kt) is about 6.5 times big- 
ger than Super-K (22.5 kt), which increases correspondingly the potential for the various 
physics applications, such as detection of relic or galactic super nova neutrinos, see Ref. [25]. 
Here we explore the possibility of a precision measurement of the solar neutrino oscillation 
parameters Am^ and 6 12 by studying the oscillations of electron anti-neutrinos emitted 
by the nuclear reactors located in the "neighborhood" of the Gd-doped MEMPHYS detec- 
tor. We will compare, in particular, the precision on Am^ and 9 12 which can be reached 
with the MEMPHYS-Gd detector with that obtainable with the Gd-doped Super-K detector 
(SK-Gd). The latter has been studied in detail in Ref. [23]. 

The water Cerenkov detectors typically do not have very good energy resolution, which 
is compensated to certain extent by their large mass. In what regards the energy resolution, 
the scintillator detectors such as KamLAND, perform significantly better. The LENA (Low 
Energy Neutrino Astronomy) detector [38] is a project for a large (~ 50 kt) liquid scintillator 
detector, to be used for studies of relic and galactic super nova neutrinos, solar neutrinos, geo- 
neutrinos, or proton decay. Since neutrinos from nuclear reactors constitute a background for 
the indicated measurements, some of the considered sites for LENA are rather far away from 
high concentrations of nuclear power plants. We consider in the following the possibility to 
place a LENA type detector in the Frejus laboratory (with many reactors relatively close by) 
and to use it for a high precision measurement of the oscillations of reactor anti-neutrinos. 

2 Reactor Neutrino Measurements with MEMPHYS-Gd and LENA 
Detectors at Frejus 



To calculate the flux of anti-neutrinos from reactors at a given position on the Earth, 
public available information on the nuclear power plants can be used [39]. A list of, and 
relevant data on reactors compiled from such sources has been kindly provided to us [40] for 
this study. To compare the reactor neutrino fluxes at Frejus with those at Kamioka, Japan, 
we consider first the effective reactor power at the detector, which is directly related to the 
total reactor v e flux reaching the detector: 

i 1 

where W* h is the thermal power of the i'th reactor, L { is the distance to the detector, and 
the the sum runs over all contributing reactors. We find that for Kamioka wj@ ani ° k!i « 
3.0MWkm" 2 /(47r), whereas at Frejus the reactor V e flux is slightly higher with H 7 e ^ ejus 
3.4MWkm- 2 /(47r). 

The average distance traveled by reactor anti-neutrinos, (L) = Wf h / Li) / (47rW e g) , 
is 188 km for Kamioka and 299 km for Frejus. Note that for Am^ = 8 x 10~ 5 eV 2 and a 
neutrino energy E u ~ 4 MeV (corresponding to the maximum of the event rate in the absence 
of oscillations, see, e.g., the first article quoted in Ref. [28]) the first oscillation minimum 
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Figure 1: Relative contribution of different reactors to the total reactor neutrino flux at Frejus and 
Kamioka as a function of the distance to the reactor. Also shown are the v & survival probabilities for 
E u = 4 MeV and 5 MeV in arbitrary units. 



of the reactor v e survival probability is at approximately 60 km. Therefore, the (L) for 
Frejus seems to be rather large for an optimal measurement of the oscillation parameters. 
However, the average distance can be misleading, and one should look at the L-distribution 
of the reactor v e flux. In Fig. Q we show the relative contribution of different reactors to the 
total reactor u e flux at Frejus and Kamioka as a function of the baseline. It turns out that 
67% of the total flux at Frejus originates from four reactors along the Rhone river located 
within a distance of 160 km from Frejus: Bugey at 115 km (25%), Saint Alban at 133 km 
(13%), Cruas at 142 km (16%), and Tricastin at 160 km (13%) 2 . Approximately 31% of the 
total flux comes from reactors distributed between 300 km and 1000 km. In our analysis we 
include 56 reactors located at a distance L < 1000 km, while the contributions of reactors at 
L > 1000 km from all around the world are summed to one "effective reactor" at 2500 km 
giving 2% of the total reactor v e flux at Frejus. 

The comparison of the L-distribution with the v e survival probability in Fig. ^ shows 
that the 4 reactors providing the dominant part of the 9 e flux at Frejus are located at 
distances which permit a rather precise study of reactor u e oscillations. For a u e energy 
E v ~ (4 — 5) MeV they are located between the first and the second survival probability 
minima, and hence spectral information should provide a powerful tool to measure the 
oscillation parameters. In the case of Kamioka the L-distribution is rather centered around 
(L) « 190 km. In Fig. the important contribution to the u e flux from the Kashiwazaki 
reactor complex located at approximately 160 km from Kamioka is clearly visible. For 
E v ~ 5 MeV this distance corresponds to the first v e survival probability maximum (see 
Ref. [27] for a detailed discussion). 

2 To use these four reactors for a measurement of B\i and Amfi has been considered previously in Ref. [29] . 



5 



We calculate the observed prompt energy spectrum by 
dN . 1 



dE p 



U Y.^J2 I dE v a{E v )UK)Pee{U,E v )R{E^E p ), (6) 



where the sum runs over the different reactors, cr(E u ) is the cross section of the detection 
reaction u e + p — > e + + n, P ee is the v e survival probability, and R(E£, E p ) is the resolution 
function relating the "true prompt energy" E^ to the prompt energy E p observed in the 
detector, where E p is determined by the initial neutrino energy, E p = E u — (m n — m p )+m e = 
E v — 0.8 MeV. We work with the total prompt energy visible in a scintillator detector 
(also if the actual detector considered is water Cerenkov) for the sake of comparison with 
KamLAND. For R(E^, E p ) we use a Gaussian resolution function with mean E^, and a 



width of 44%/ ^E p [MeV] for MEMPHYS-Gd/SK-Gd and 10%/y/E p [MeV] for LENA [38]. 
The energy resolution for MEMPHYS-Gd and SK-Gd is similar to the one reported by 
Super-K for the solar neutrino analysis (see Fig. 15 of the second paper in Ref. [2]). 

In Eq. ©, 4>i(E u ) denotes the flux of v e emitted by reactor i, which is given by 

<t>i(Ev) = W?J2wMEv), (7) 

where i = 235 U, 238 U, 239 Pu, 241 Pu, labels the most important isotopes contributing to the 
P e flux, fi is the relative contribution of the isotope I to the total reactor power, and Eg is 
the energy release per fission for the isotope I. In Eq. (JZj), (fie(E u ) is the (energy differential) 
number of neutrinos emitted per fission by the isotope £, and we adopt the parameterization 
for the <pg(E u ) from Ref. [41]. For the fe we take a typical isotope composition in a nuclear 
reactor of [5] 235 U : 238 U : 239 Pu : 241 Pu = 0.568 : 0.297 : 0.078 : 0.057, and we assume these 
ratios to be equal for all reactors. 

In the calculations we use the 3-neutrino oscillation survival probability P ee which de- 
pends, in particular, on sin 2 9 is (see, e.g., [10]), and take into account the (small) Earth 
matter effect. In the case of absence of oscillations, P ee = 1, the number of events above 
a threshold E thr is given by iV noosc = AfW c sC, where W e s has been defined in Eq. (J3J) and 
C is an integral depending only on E^ t . For MEMPHYS-Gd we use a threshold for the 
prompt energy E^ T = 3.0 MeV (which corresponds to the value of 2.5 MeV for the positron 
energy given in Ref. [25]), whereas for LENA we use E^ = 2.6 MeV to eliminate the back- 
ground from geo-neutrinos, as in the KamLAND oscillation analysis [5,6]. To determine 
the normalization constant M in Eq. © we use the prediction for N QOOSC in KamLAND [5], 
and then we scale it for each experiment taking into account that M is proportional to the 
measurement time and the number of free protons in the detector, as well as the different 
values of W e s and E t hr- In Tab. Q]we summarize the most important characteristics of the 
considered detectors as simulated in our analysis, and we give the expected number of events 
in case of no oscillations. 

To test the sensitivity of the experiments we divide the prompt energy spectrum in 
Eq. © into 20 bins between 3 MeV and 12 MeV for MEMPHYS-Gd and SK-Gd, and into 
25 bins between 2.6 MeV and 10 MeV for LENA 3 . The data is simulated using as "true 

3 For low statistics data samples such as the present KamLAND one, a likelihood analysis [42] or equal 
bins in l/E p [21] allows to extract an optimum of information. In the cases under study, however, the 
numbers of events are relatively high. Therefore, we simply take equal bins in E p , with a bin size sufficiently 
smaller than the energy resolution to make sure that no information is lost due to the binning. 
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experiment 


na. mass 


free protons 


771 

-^thr 


event s/yr 


enerj 


$y resol. 


MrL/Mr ri Y o-<ord 


147. U kt 


n o x , i n33 

9.8 x 10 


3.0 MeV 


59 980 


A A 0/ 1 

44%/ v 


'hp [MeV J 


LENA 


44.0 kt 


2.3 x 10 33 


2 6 MeV 


16 670 


10%/a 


/ E fMe VI 


SK-Gd 


22.5 kt 


1.5 x 10 33 


3.0 MeV 


8 000 


44%/ v 


IE V [MeV] 


KamLAND 


0.41 kt 


3.5 x 10 31 


2.6 MeV 


216 


7.5%/- 


sjE v [MeV] 



Table 1 : Summary of the input characteristics of the detectors MEMPHYS-Gd, LENA, and SK-Gd, used 
in our analysis. For comparison we show also the corresponding values for KamLAND [5]. The number of 
events/yr is calculated for no oscillations and using the reactor flux at Frejus for MEMPHYS-Gd and LENA, 
and at Kamioka for SK-Gd. 



I 1 I 1 




J i I i I i I i I i I i i 

3 4 5 6 7 8 9 10 11 12 
E [MeV] 



Figure 2: The ratio of the event spectra in positron energy in the case of oscillations with Amfi = 
7.9 x 10~ 5 eV 2 and sin 2 6*12 = 0.30 and in the absence of oscillations, determined using one year data of 
MEMPHYS-Gd and LENA located at Frejus. The error bars correspond to ler statistical error. 

values" Am^ = 7.9 x 10~ 5 eV 2 and sin 2 9 12 = 0.30. The latter correspond to the present 
best fit point obtained in a global neutrino oscillation analysis [21]. Then a x 2 -analysis is 
performed to determine the allowed regions and the precision with which these parameters 
can be determined from the simulated data. In Fig. El we show the ratio of the number of 
events calculated by taking into account v e oscillations with parameters indicated above to 
the number of events in the absence of oscillations, binned in prompt energy. The error 
bars correspond to la statistical errors for one year of (simulated) MEMPHYS-Gd and 
LENA data. The large number of events leads to a very precise measurement of the energy 
spectrum. The oscillatory signal in the spectrum is less pronounced in the MEMPHYS-Gd 
spectral "data" than in the analogous LENA "data" due to the worse energy resolution of 
the water Cerenkov detector. Nevertheless, as a consequence of the relatively high statistics, 
a clear signal of spectral distortion can still be observed with MEMPHYS-Gd. In the case of 
LENA spectral "data", an event maximum is clearly visible at E p = (3.5 — 4.0) MeV, which 
originates from the first oscillation maximum of the survival probability at L = 160 km (see 
Fig. Q). In both cases the spectral information is crucial for the precise determination of the 
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systematic value 

overall normalization (fully correlated) 5% 

thermal power of each reactor (uncorrelated) 2% 

energy scale uncertainty 0.5% 

prompt energy spectrum tilt 2% 

reactor neutrino flux Ref. [41] 



Table 2: Systematical uncertainties and the default values adopted in our analysis. 

oscillation parameters. 

In the statistical analysis we take into account various systematical uncertainties as listed 
in Tab. |21 We include a 5% error on the overall normalization (e.g., from the uncertainty 
on the fiducial mass), a 2% uncertainty on the thermal power of each reactor (uncorrelated 
between the reactors), and the uncertainty on the anti-neutrino spectra <j>e(E v ) (normaliza- 
tion and shape) according to Ref. [41]. We take into account an uncertainty of 0.5% in the 
energy scale calibration of the detector. This value is motivated by the numbers given for 
the Super-K solar neutrino analysis (see second reference in [2]) and for the Double-Chooz 
reactor experiment [43]. Furthermore, we include a linear tilt in the prompt energy spectrum 
of 2%, i.e., we allow the event number in the highest energy bin to shift by 2% with respect 
to the event number in the lowest energy bin with a linear interpolation for the intermediate 
bins. In the following section we will discuss in some detail how much our numerical results 
depend on the values adopted for the systematic uncertainties. In addition to these system- 
atics we have tested also the effect of an uncertainty on the isotope compositions fi defined 
in Eq. (JZJ). We have performed an analysis allowing the fg to vary independent for each 
reactor within 5%, and found that the impact on the sensitivity to the neutrino oscillation 
parameters is negligible. Therefore, we keep the fe fixed in our standard analysis which 
significantly reduces the calculation time. 

3 Precision of the Determination of Neutrino Oscillation Para- 
meters 

In this Section we present results on the precision which can be reached in the measure- 
ment of Amlx and sin 2 9^ in the MEMPHYS-Gd and LENA experiments. Our findings are 
summarized in Fig. El and Tab. 01 where we compare the results which can be obtained using 
one year of data from MEMPHYS-Gd, LENA, and SK-Gd, with the present constraints from 
the global solar and KamLAND data [19,21]. In the table we give the relative uncertainty 
at 3a defined as 

^.upper flower 

spread(x) = ; , (8) 

1 v ' flipper _|_ flower ' *> ' 

where a; u PP er ( a ; lower ) j s the upper (lower) bound of the quantity x at 3a. 

We find that the water Cerenkov detector MEMPHYS-Gd and the scintillator detector 
LENA can provide very similar high precision determinations of Am^ and sin 2 9i2- The 
lower mass of LENA is compensated by its better energy resolution. Already with one year 
of data, an uncertainty smaller than 3% at 3a can be obtained in the determination of 
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spread(Am2 1 ) spread(sin 12 ) 



time 


1 yr 


7 yr 


i yr 


7 yr 


SK-Gd 


6.0% 


2.8% 


36.6% 


18.6% 


MEMPHYS-Gd 


2.9% 


1.4% 


20.0% 


13.2% 


LENA 


2.5% 


1.2% 


18.0% 


9.8% 


solar + KamLAND 




11.3% 


24. 


9% 



Table 3: The 3er uncertainty as defined in Eq. (JSJ in the determination of Am^ and sin 2 812 in the 
experiments MEMPHYS-Gd and LENA at Frejus and SK-Gd at Kamioka after 1 and 7 years of data taking. 
For comparison we show also the current uncertainties in Anijj and sin 2 #12 [21]. 




0.1 0.2 0.3 0.4 0.5 10 20 30 40 
sin 2 A% 2 



Figure 3: The accuracy of the determination of Am 2 ^ and sin 2 612, which can be obtained using one year 
of data from MEMPHYS-Gd and LENA at Frejus, and from SK-Gd at Kamioka, compared to the current 
precision from solar neutrino and KamLAND data. We show the allowed regions at 3<r (2 d.o.f.) in the 
Am|i — sin 2 612 plane, as well as the projections of the x 2 for each parameter. 

Amjj, while sin 2 # 12 can be determined with an error of about 20% at 3a. This precision is 
approximately by a factor two better than the precision that can be reached with one year 
of data from SK-Gd. The better precisions which can be obtained with the MEMPHYS-Gd 
detector compared to those that can be obtained with the SK-Gd detector are a consequence 
of the larger fiducial mass of MEMPHYS-Gd. As follows from Tab. El and Fig. 1 year 
of data from MEMPHYS-Gd and 7 years of data from SK-Gd (yielding approximately the 
same numbers of events in the two detectors, see Tab. allow a determination of Am^ and 
sin 2 #12 with similar precisions. This shows, in particular, that both locations, Frejus and 
Kamioka, are very similar in what regards the power of the "surrounding" reactors and their 
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0.1 0.2 0.3 0.4 0.5 10 20 30 40 
sin 2 1? A% 2 



Figure 4: The allowed regions in the Am|i — sin 2 612 plane, obtained at 3a (2 d.o.f.) from 147 kt yr 
and 147 x 7 kt yr of MEMPHYS-Gd, and from 22.5 kt yr and 22.5 x 7 = 157.5 kt yr of SK-Gd data. The 
projections of the \ 2 f° r each parameter are also shown. 

distance distribution. Thus, the two locations are equally suitable for the high precision 
measurements under discussion. 

Ultimately, 7 years of data from MEMPHYS-Gd (LENA) would allow a determination 
of Amji and sin 2 612 with uncertainties of approximately 1.4% (1.2%) and 13% (10%) at 
3a, respectively. This precision is comparable to the precision which can be reached in 
the determination of the atmospheric neutrino oscillation parameters Am^ and sin 2 #23 by 
studying disappearance in the superbeam experiments T2HK in Japan or SPL from CERN 
to MEMPHYS (see e.g., Ref. [37] for a recent analysis). Hence, the reactor measurement 
could complete the program of the high precision determination of the parameters responsible 
for the leading solar and atmospheric neutrino oscillations. 

We have also investigated whether the uncertainty in the knowledge of #13 can have any 
effect on the precision of Am 2 ^ and Qyi determination in the experiments MEMPHYS-Gd and 
LENA under discussion. We show in Fig. El the results of a three- flavour neutrino oscillation 
analysis of 1 year simulated data in MEMPHYS-Gd, assuming that the true value of #13 is 
zero. In the analysis sin 2 #13 was allowed to vary freely, however the information available 
at present has been included in the fit by adding the X 2 ($i3) obtained from current global 
data in Ref. [21]. In the panel for the sin 2 #12 — Am^ projection we show also the allowed 
regions, obtained in a two-flavour analysis, i.e., for # 13 = 0, with black contours. Indeed, 
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"7 I i i i I i i i I i i i I i i i I I i I i I i I i 

0.1 0.2 0.3 0.4 0.5 0.05 0.1 0.15 0.2 

2 2 

sin 9 12 sin 20 13 

Figure 5: Allowed regions from a three-flavour neutrino oscillation analysis of 147 kt yr data in 
MEMPHYS-Gd at 90%, 95%, 99%, and 99.73% CL (2 d.o.f.), projected onto the three different 2-dimensional 
parameter planes. The data are simulated for #13 = 0. In the analysis #13 was varied freely taking into ac- 
count the constraint from the current global data as obtained in Ref. [21]. The allowed regions in the 
sin 2 #12 — Am 2 ,! plane from a two-flavour analysis with #13 fixed to zero (indicated with black contours) are 
identical to those obtained in the three-flavour analysis. 

the regions within the black contours are practically identical to the shaded/colored regions 
corresponding to the three-flavour case. Therefore we can conclude that the uncertainty 
in the knowledge of (9 13 does not affect the Am^ and # 12 measurements. As is visible in 
Fig. El there are no correlations of the leading parameters with #13, since Am 2 ^ and #12 are 
determined primarily from the data on the shape of the spectrum, which does not depend 
on 6*13. We have checked that the above conclusions concerning the 6*i 3 -effects hold also for 
LENA and SK-Gd detectors. 

Let us note that the sensitivity of MEMPHYS-Gd to #13 on its own is rather poor. The 
constraint on sin 2 9 13 appearing in Fig. El corresponds just to the bound from present global 
data, which is included in the analysis. Hence, the sensitivity of MEMPHYS-Gd is worse 
than the present bound. This is a consequence of the fact that a non-zero #13 basically leads 
to a rather small (energy-independent) suppression of the total v e flux, which is unobservable 
due to the relatively large uncertainties in the overall normalization. 

Since we are discussing here very high statistics measurements, a careful investigation 
of the impact of systematical uncertainties is necessary. In Fig. El we show how our results 
for the accuracies of the oscillation parameters depend on the assumptions adopted for the 
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MEMPHYS-Gd (lyr) ^■LENA(lyr) E^S SK-Gd (7yr) 



All systematics 
Reactor power 
Energy scale 
Spectrum tilt 
Normalization 

0.6 0.8 1 1.2 1.4 0.4 0.7 1 1.3 1.6 

2 2 

impact on Am J2 impact on sin 

Figure 6: Impact of systematical uncertainties. We show the relative change in the 3<r intervals of Am^ 
and sin 2 612 obtained by switching off (left edges of the bars) and doubling (right edges) all systematics 
simultaneously, the uncertainties on the thermal power of the 4 most important reactors ("Reactor power"), 
the energy scale uncertainty, the prompt energy spectrum tilt, and the overall normalization error. 




systematical errors. In particular, we consider the impact of uncertainties on the thermal 
power of the 4 most important reactors, the energy scale uncertainty, the prompt energy 
spectrum tilt, the overall normalization error, as well as all systematical errors in total. To 
check the impact of these uncertainties we show in Fig. El the ratios S /5 st d and 5 2 /<W, where 
So(2) is the 3<T range for Am| or sin 2 9\2 if the systematical error of interest is set to zero 
(is doubled with respect to its standard value), and 5 s t& is the 3a range using our standard 
values according to Tab. 121 

The left edges of the bars in the row denoted by "All systematics" in Fig. El correspond 
to statistical errors only. In this ideal case the Am| x accuracy is improved by about 20% for 
MEMPHYS-Gd, 10% for LENA, and 30% for SK-Gd with respect to our standard choice for 
the systematics, whereas the precision on sin 2 9 12 is improved for all experiments by more 
than a factor of 2. For the Am| measurement the individual systematics have only a minor 
impact (with the exception of a ~ 20% effect of the energy scale uncertainty in SK-Gd). 
For the measurement of sin 2 9 12 the overall normalization and the energy scale (especially 
for MEMPHYS-Gd) are important. Note that the uncertainty on the thermal reactor power 
has a negligible impact on the accuracies. Therefore, it seems not to be possible to improve 
the precision on Am^ and 9 12 by installing near detectors close to the reactors dominating 
the v e flux. 

In summary, systematic uncertainties are an important factor in the experiments under 
consideration. Especially the determination of the mixing angle depends on the values 
of systematic errors. The overall effect emerges from an interplay of the various sources 
of uncertainties included in our analysis, and therefore, to obtain an improvement in the 
precision of the oscillation parameters several of the systematic errors listed in Tab. Elshould 
be decreased. 
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4 Conclusions 



We have investigated the possibility of a high precision determination of the solar neutrino 
oscillation parameters Am^ and 9 12 in a long-baseline reactor neutrino experiment, located 
in the Frejus underground laboratory. Approximately 67% of the total reactor V e flux at 
Frejus originates from four nuclear power plants in the Rhone valley, located at distances 
between 115 km and 160 km from Frejus. The indicated baselines are particularly suitable 
for the study of the u e oscillations driven by Ara^ — they are similar to those exploited in the 
KamLAND experiment in Japan. Approximately 31% of the total flux v e at Frejus comes 
from reactors distributed between 300 km and 1000 km from Frejus. In our analysis we 
include 56 reactors located at a distance L < 1000 km, while the contributions of reactors at 
L > 1000 km from all around the world are summed to one "effective reactor" at 2500 km 
giving 2% of the total reactor u e flux at Frejus. The Frejus underground laboratory is under 
consideration as a possible site for a mega ton scale water Cerenkov detector MEMPHYS 
which, among other physics applications, may serve as a far detector for a neutrino beam 
produced at CERN. In the present article we have assumed that the water of one module of 
MEMPHYS having a fiducial mass of 147 kt, is doped with 0.1% Gadolinium (MEMPHYS- 
Gd), which will allow, in principle, a high precision study of reactor V e oscillations. As an 
alternative detector technology, we have considered a 50 kt scale liquid scintillator detector, 
as discussed in the LENA proposal, which can be viewed as a considerably larger version of 
the present KamLAND or Borexino detectors. 

The analysis performed by us shows that each of the two detectors — MEMPHYS-Gd and 
LENA, if placed at Frejus, would allow a very precise determination of the solar neutrino 
oscillation parameters Am^ and sin 2 #12: with one year of reactor V e data taken at Frejus 
(by any of the two detectors), the 3a uncertainties on Am^ and sin 2 9 12 can be reduced 
respectively to less than 3% and to approximately 20%. The Gadolinium doped Super- 
Kamiokande detector (SK-Gd) in Japan can reach a similar precision if the SK/MEMPHYS 
fiducial mass ratio of 1 to 7 is compensated by a longer SK-Gd data taking time. Several 
years of reactor V e data collected by MEMPHYS-Gd or LENA would allow a determination 
of Am^ and sin 2 6*12 with uncertainties of approximately 1% and 10% at 3<r, respectively. 
We have shown also that the uncertainty associated with the CHOOZ mixing angle 13 has 
practically no impact on the measurements of the solar neutrino oscillation parameters in the 
experiments discussed by us, and we have investigated in some detail the effects of various 
systematical uncertainties on the precision of the determination of Am^ and sin 2 9\2 in these 
experiments. 

The accuracies on the solar oscillation parameters, which can be obtained in the high 
statistics experiments considered here are comparable to those that can be reached for the 
atmospheric neutrino oscillation parameters Am^ and sin 2 #23 in future long-baseline su- 
perbeam experiments like T2HK in Japan or SPL from CERN to MEMPHYS. Hence, such 
reactor measurements would complete the program of the high precision determination of 
the leading neutrino oscillation parameters. 
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